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ABSTRACT 

We present the first spatially resolved polarized scattered light H-band detection of the DoAr 28 
transitional disk. Our two epochs of imagery detect the scattered light disk from onr effective inner 
working angle of 0710 (13 AU) out to Cf.^50 (65 AU). This inner working angle is interior to the location 
of the system’s gap inferred by previous studies using SED modeling (15 AU). We detected a candidate 
point source companion 1708 northwest of the system; however, our second epoch of imagery strongly 
suggests that this object is a background star. We constructed a grid of Monte Carlo Radiative Transfer 
models of the system, and our best fit models utilize a modestly inclined (50°), 0.01 Mq disk that has 
a partially depleted inner gap from the dust sublimation radius out to ~8 AU. Subtracting this best fit, 
axi-symmetric model from our polarized intensity data reveals evidence for two small asymmetries in the 
disk, which could be attributable to variety of mechanisms. 
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1. INTRODUCTION 

Detailed observations of the spectral energy distribu¬ 
tions (SEDs) of young stellar objects (YSOs) have revealed 
numerous systems with deficits of near-infrared flux com¬ 
pared to primordial YSOs. These pre-transitional and 


gaps and holes i 

n their disks (Strom et al. 1989 Espail- 

lat et al.||2007a|b 

). Both infrared scattered light imagery 

dFukagawa et al. 

2006 Thalmann et al. 2010 

Hashimoto 

et al. 120111 Canovas et al.|2013 Avenhaus et a 

.120141) and 

sub-millimeter observations d Andrews et al. 2011 Brud- 

erer et al.|2014 a 

^an der Marei et al. 112014) have contirmed 


disk systems, and revealed additional sub-structure that 
could provide information about the mechanism responsi- 


ble for clearing regions in these disks (van der Marel et al. 

M4|. 


Numerous mechanisms have been suggested to explain 
the clearing and partial clearing associated with pre- 
transitio nal and transitional disk syste ms, including grain- 
growth (Dullemond fc Dominik 20051, photo-evaporation 
(Clarke et al. 2001[ ), disk instability ( Papaloizou 2007 1, 
and perturbation from planetary companions (Zhu et al. 


2011[). These mechanisms can also influence other regions 


of these disks such as the inner walls and surface struc¬ 
tures. Fo r example, grain grow th can create a rounded 
gap edge (Birnstiel et al. 20121 whereas planets can in¬ 
flate gap edges, creating a wall that can shadow parts of 


the outer disk or induce signif icant back scattering (Jang- 
Condell fc Turne'r]|2012 2013). Grain growth should also 


induce very bright sub-millimeter features in systems ex¬ 
hibiting gaps in the near-infrar ed, although this is not seen 
in Oph IRS 48 or SAG 206462 ( Bruderer et al.|2014 Perez 


et al. 2014). Additional morphological features such as spi- 
ral arms have been detected in these disks, and could be 
attributable to stellar and sub-stellar co mpanions or grav¬ 
itational instabilities (Bate et al. 2003). Planet-induced 
perturbations should co-rotate with the planet; hence, de¬ 
tecting the rotation of these structure could help to dis¬ 


tinguish the origin of some d isk structures (Hashimoto et 
al. 2011 Lomax et al. 2015). Additional morphological 


futures observed in transitional disks, such as dust traps 


(Bruderer et al. 2014), can also be used to constrain the 
mechanism responsible for sculpting the spatial distribu¬ 
tion of gas and dust in transitional disks. Clearly, a first 
step that is needed to assess the clearing mechanism is to 
fully constrain the spatial distribution of small and large 
dust grains, as well as the gas, in individual disk systems. 


ous transitional disks (jHashimoto et aL|2012 

Muto et al. 

2012[ Grady et al. 20131 Follette et al. 201S 

Kusakabe 

et al. 1120121 IHashimoto et al.pOll Thalmann et al.||2010 

Tsukagoshi et al.l|2014 

Takami et al.||2013 Canovas et al. 

2013 Avenhaus et al. 

2014 Hashimoto et al.|2015) have 


eluding large programs like the Strategic Exploration of 
Exoplanets and Disks with Suba ru (SEEDS) high-contrast 
imaging survey ( Tamura~||2009 ). These observations have 


provided scattered light conhrmation of the gapped nature 


of transitional disks (Thalmann et al.]|2010), revealed the 


presence of spiral structures (Muto et al. 2012 Grady etal. 


2013 

Currie et al.|2014|, detected likely non-axisymmetric 

inner 

disks (Kusakabe et al. 

2012 

Takami et al. 

2013 

1, 


found differences m some cases about the distribution of 


al. 

2012 

Dong et al. 2012| 

, and the presence of a planet 

within a gap (ICurrie et al. 

2014 

I. 


DoAr 28 is a K5-type object in the p Ophiuchi associ - 
ation located at a distance of ^139 pc (Mamajek 2008), 
that has been identified as a transitional dis k based on its 
SED (McClure et al. 2010 Kim et al.||2013 ). The system 
is ob served to be aHive ly accreting at a rate of 4 * 10“® 


dKeane et al. 

year ^ 


2014), and has a disk with an inferred 
15 AU (Kim et al.|2013 ) from analysis 


outer gap radius of 
of its SED. 

In this paper, we present the first scattered light detec¬ 
tion of the DoAr 28 transitional disk, in the H-band. In 
Section 2, we will discuss the scattered light observations 
and data reduction. In Section 3, we present our analysis 
of the resolved disk, constrain the presence of co-moving 
point sources, and model these data using Monte Carlo 
models. Finally, we discuss the implications of our results 
in the broader context of other resolved transitional disk 
systems in Section 4. 

2. OBSERVATIONS AND REDUCTIONS 

DoAr 28 was observed in two epochs on 2012 July 9 
(2012 epoch) and 2014 June 9 (20 14 epoch) as part of 


the SEEDS survey using HiCIAO (Hodapp et al. 


2008) 


in the H-band. The data were obtained in quad Polar 
ized Differential Imaging (qPDI) mode at four wave-plate 
positions (0°, 22.5°, 45°, 67.5°) in Angular Differential 
Imaging (ADI) mode. Each observational frame contains 
four sub-images with each sub-image having a field of view 
of 570 by 570, with a pixel scale of 9.5 mas pixel“^ and 
a FWHM of (y.'136 for the 2012 epoch and (X'lOl for the 
2014 epoch. The total ADI field rotation achieved was 
19.5° for the 2012 epoch data and was 21.3° for the 2014 
epoch data. 

2.1. HiCIAO Data Reduction 

The data reduction process we employed for extracting 
polarized intensity (PI) disk images utilized the double 


et al. 

(2011 

each 

rame c 


(2011). To briefly review, the four sub-images of 


images, which can be summed and subtracted from their 
90° counterparts to create -Q, -l-Q, -U, and -|-U images. 
The Q and U frames were then rotated into a common 
orientation, corrected for instrumental polarization, and 
summed to create final Q and U images. Note that as 
our 2012 epoch imagery were obtained under non-optimal 
conditions, we only utilized the best 64 of the 76 observed 
frames for our summed imagery. The final PI images are 
computed using PI = + U^. As previously noted by 


Hashimoto et al. (2012), the PSF convolved by seeing is 
not perfectly corrected by the AO-188 system, which pro¬ 
duces a residual polarized halo. We computed an artificial 
halo following the procedure outlined in [Hashimoto et ah] 
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(2012), scaled this to the observed aperture polarization 
of the system (P = 0.869% ± 0.012%), and subtracted it 
from the PI image to create final PI imagery. 

We searched for point-source companions to DoAr 28 


using the ACORNS-ADI software package (Brandt et al. 


2013). We treated each of the four sub-images noted above 


as its own angular differential imaging sequence (ADI, 
Marois et al. 2006) and reduced them using the LOCI 
algorithm (Lafreniere et al. 2007) with the standard reduc¬ 
tion parameters from Brandt et al. (2013). This yielded 
four residual images, each corrected for partial flux sub¬ 
traction. We then averaged these PSF-subtracted images 
to produce a single high-contrast image. We computed 
the standard deviation in annuli on this combined image 
to produce a contrast map and searched for 5cr compan¬ 
ions. We found one companion candidate as described in 
Section 3.2; follow-up data showed it to be an unrelated 
background star. 


2.2. SMA Observations and Reduction 

We observed DoAr 28 with the Submillimeter Array 
(SMA) on 2011 March 16, using the Compact Configu¬ 
ration with six of the 6 m diameter antennas at 230 GHz 
(1.3 mm) with a full correlator bandwidth of 2 GHz, for 
a total integration time of 63 minutes. Calibration of the 
visibility phases and amplitudes was achieved with obser¬ 
vations of the quasar 3C 279, at intervals of about 20 min¬ 
utes. Observations of Titan provided the absolute scale 
for the flux density calibration. The data were calibrated 
using the MIR software package]^ We detected DoAr 28 
with a flux density of 68.6±1.9 mJy. The double sideband 
system temperatures were 110 to 170 K. 


2.3. APO Observations and Reduction 


where the radial profile measured along the disk major axis 
(Figure]^ exhibited a clear deviation from the disk domi¬ 
nated scattered light power law behavior seen in the outer 
disk. Using this criterion, we determined the inner working 
angle to be 0'.'17 (22 AU) for the 2012 epoch data and Ol'lO 
(13 AU) for the 2014 epoch (FigurejS. The power law that 
characterizes the radial profile distribution of the scattered 
light flux along the disk major axis extends out to ^U.'SO 
(65 AU), defining the outer edge of the scattered light disk, 
before becoming clearly dominated by noise. The detected 
disk appears to be continuous, with no significant gaps or 
holes clearly visible. The two epochs of imagery appear 
similar to each other, modulo potential differences aris¬ 
ing from the factor of 1.4 better FWHM achieved in the 
2014 epoch imagery. For the imagery with a slightly larger 
FWHM, it is conceivable that a small amount of disk flux 
could be spread out into the PSF halo, and be subtracted 
out in the process described in Section |2.1[ Nevertheless, 
the overall surface brightness of the disk is the same at 
both epochs and the radial surface brightness power law 
measured along the major axis for both epochs is similar 
(-2.47 for 2012; -1.84 for 2014). We do note the poten¬ 
tial presence of a slight curl in the SW region of the 2014 
epoch scattered light disk, and further discuss this feature 
in Section |4] 

We determine the inclination of the disk to be ^50° 
(Table [^. This was found by subtracting a suite of disk 
models Tor a range of inclination angles, described in Sec¬ 
tion and identifying the inclination which yielded the 
smallest residuals (Figure [^. Note that this residual im¬ 
age (Figure]^ exhibits clear deviations from axisymmetry, 
with a deficit of scattered light present along the north¬ 
ern side of the major axis. We will be further discuss this 
non-axisymmetric structure in Section]^ 


DoAr 28 was observed using the ARC Echelle Spectro¬ 
graph (ARCES; [Wang et al. 2003) at the Apache Point 
Observatory (APO) 3.5m telescope on 2014 June 19, yield¬ 
ing a R~31,500 spectrum covering the spectral range of 
^3,600-10,OOOA . The data were reduced using standard 
IRAF techniques. We extracted the order containing Ha 
and continuum normalized these data, enabling us to char¬ 
acterize the line strength discussed in Section [T3| 


3. ANALYSIS 

3.1. Scattered Light Imagery 

Our two epochs of scattered light PI imagery of DoAr 28 
are shown in Figure Significant scattered light around 
the central star is clearly seen in both epochs. The direc¬ 
tion and relative intensity of the polarization vectors de¬ 
rived from these data exhibit a clear centrosymmetric be¬ 
havior around the central star (Figure]^, confirming that 
this signal arises from scattering off of circumstellar mate¬ 
rial. These data exhibit less evidence of centro-symmetry 
about the minor axis, which is much less resolved than 
the major axis, suggesting that the level of our residual 
polarized halo correction might be incomplete. Our pri¬ 
mary analysis of these data will focus on the measured 
polarized intensity along the major axis of the system. 

Since these data were not observed with a coronagraph, 
we defined the effective inner working angle as the location 

http://www.cfa.harvard.edu/~cqi/mircook.html 


3.2. Point Source Detections 

We identified a candidate point source companion 1/09 
northwest of DoAr 28, with an i7-band contrast of 9.5 
magnitudes, in our July 2012 data. Our second epoch im¬ 
ages from June 2014 indicate that this object is almost 
certainly a background star. DoAr 28 was too faint for a 
proper motion from the Tycho satellite, but is part of to 
the same star forming region as p Oph. We assume that it 
shares p Oph’s proper motion of (—5.5 ± 0 .9, —21.7 ± 0.9) 
mas/yr in RA and Dec ( van Leeuwen|2007 ). At a distance 


of ~^140 pc, 1 mas/yr corresponds to ~0.7 kms similar 


TW Hva (Perrvman et al. 

1998 

Jones 

1970 

Makarov & 

Fabricius|2001 ). and somewhat less than the ~2- 

3 kms"^ 

of Orion (jJones & Waiker|1988 

Ffiresz et al. 2008 

I. Figure 


background track assuming common proper motion with 
p Oph. If the point source were co-located with DoAr 28, 
their relative velocity would be at least 17 kms“^, much 
too high for the system to be bound at its observed loca¬ 
tion. 

Figure shows our sensitivity limits, assuming a dis¬ 
tance of ~^140 pc and converting contrast to absolute mag¬ 
nitude, and neglecting extinction (estima ted to be ~^0.4 
magn itudes at H-band, using A(v) = 2.3 ([McClure et al. 


2010), Rt, = 3.1, and the Cardelli, Clayton, &: Mathis 
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(1989) extinction relations), making these sensitivity es¬ 


timates slightly optimistic. It is more difficult to inter¬ 
pret these results as mass limits, given DoAr 28’s extreme 


young planets dSpiegel & Burrows|20I2 

Allard et al. 12011 

Marley et al.||2007| 

I'ortney et al.||2008 

Our yEEDS ob- 

servations reach a . 

imiting 71-band absolute magnitude of 


^13.5 at a projected separation of 100 AU. At an age of 5 
Myr, this corresponds to ^4 Mjup in the BT-Settl models, 
but anywhere from ~4 Mjup up to the deuterium- burning 
limit of ^13 Mjup in the Spiegel & Burrows (2012) models 


(SB12), depending on the initial entropy. We note that the 
high-mass end of this range requires a very cold start. As¬ 
suming an initial entropy midway between the minimum 
and maximum values of the SB12 models, the SB12 mass 
limits are only slightly higher than the BT-Settl limits. 

There are theoretical reasons arguing against the for¬ 
mation of substellar co mpanions below ^5 Mjup by direct 


gravitational collap se (Low & Lynden-Bell 1976 Bate et 
al.|2003 Bate|2009 ). Our SEEDS imagery rule out such a 
companion, which would necessarily form hot and beyond 


companion, wtiicli would necessarily form hot and beyond 
~80 AU. While a less massive core-accretion planet is un¬ 
likely to form at such a wide separation (DoAr 28’s disk 
only extends to 70 AU), our data cannot rule out a Jovian 
planet scattered into a wide or unbound orbit. 

3.3. Radiative Transfer Modeling 

We modeled the DoAr28 system with the 
HOCHUNK3D Monte Carlo Radiative Transfer (MCRT) 
code as described in Whitney et al. (|2013 |. 
HOCHUNK3D is similar to other codes ( Min et ^ 2009' 


Wolf fc Hillenbrandl |2003 |Dullemond fc Dom inik 


Pinte et al. 2006 Robitaille 2011 D’Alessio et al. 


2004 


2006) 


and has a long history of being used to constrain the dust 
distribution in protoplanetary systems. The latest version 
of HOCHUNK3D decouples the small and large dust grain 
distributions, allowing settling of dust to be incorporated. 
This implementation can be thought of as utilizing over¬ 
lapping disks with different dust grain size distributions. 
Dust density distributions i n HOC HUNK3D are adopted 
from Shakura & Sunyaev ( 1973[ ) and are characterized 
by a radial power law (a), and a vertical gaussian distri¬ 
bution (/3), as given by equation]^ where r is the radial 
component in cylindrical coordinates and 2 is the height of 
the disk from the mid-plane, h is the scale height defined 
in equation which is normalized at a defined radius of 
100 AU. Deviations in the dust density distribution, such 
as gaps, spiral arms, and warped disks, can also be fully 
parameterized with HOCHUNK3D. 


p oer “ exp 


-h 


h (xr ^ 


( 1 ) 


( 2 ) 


Whitney et al. (2013) and references therein describe 
the full radiative transfer of the HOCHUNK3D code. To 
briefly summarize, the code uses a Henyey-Greenstein 
scattering phase function and includes parameters for 
forward-scattering and albedo calculated from the adopted 
dust grain model. The dust models we used are described 
below. Temperature corrections utilized the Lucy method 
with a maximum number of six iterations (Lucy 1999). 


We utilized 5*10® photons for our broad exploration of 
MCRT parameter space, and 5*10^ photons for each of 
our runs where we compared both the observed SED and 
PI imagery against the models. 

As is true with many MCRT codes, the large number of 
free parameters exceed the number of data points leading 
to model degeneracies. Since we lacked spatially resolved 
sub-millimeter observations that trace the radial distribu¬ 
tions of large grains we assumed that the small and large 
grain disks had the same a and /3 . We included a wall 
along the disk edge, as walls are though t to be common 
in transitional disks (Calvet et al. 2005 Espaillat et al. 


2007a ). We also assumed that there was a negligible en¬ 
velope, no warping of the disk, and that the disk was az- 
imuthally symmetric. We adopted dust parameters from 
Wood et al. dust model 1 (2002) for our large dust grains, 
which are composed of amorphous carbon and silicon rang¬ 
ing in sizes up to 1 millimeter. The small dust grain model 
we ado pted w as the average galactic ISM model from Kim 
et al. (1994). We generated a grid of 270 models and 
identified a broad range of parameters consistent with the 
observed SED in Table [T] 

Our model SED is consistent with the observed pho¬ 
tometry of DoAr 28, as shown in Eigure with the new 
SMA data point helping to constrain the large dust grain 
disk. We explored a range of gap sizes and how these in¬ 
fluenced the resultant SED and images, and demonstrate 
below that a gap size of ^8 AU best represents our data. 
We also found that the accretion rate and the gap den¬ 
sity parameter had similar effects o n the SED T| We 


(2014), 


adopted the accretion rate quoted by Keane et al. 

4.0* 10“® Msl ^ which led us assume a gap density param¬ 
eter of 5*10“® (Table [^. Note however that our analysis 
of our new spectroscopic observations of DoAr 28 revealed 
a H-alpha equivalent width (30±0.lA) that was different 


Keane et al. 


(2014). This sug- 


than the 36 A reported in 
gests that the system likely exhibits a variable accretion 
rate. Thus, the upper and lower bound values of the gap 
density shown in Table [l] are more representative of the 
system. 

After using the observational SED of DoAr 28 to con¬ 
strain our MCRT model parameters, we next used the 
observed surface brightness of our PI imagery to further 
constrain these parameters. Specifically, we measured the 
surface brightness of our PI imagery using a 4-pixel wide 
aperture along the major axis of the disk, and compared 
this to the H-band PI surface brightness predicted by our 
models (Figure®. This iterative process enabled us to ar¬ 
rive at our finaladopted model parameters listed in Table 
[T] including the adoption of the gap size at ^8 AU. We 
used the model image subtracted from the PI images to 
find the inclination of 50°. To search for potential devi¬ 
ations from axi-symmetry in our data, we subtracted our 
(axi-symmetric) model, scaled to the peak intensity of the 
observed PI disk, from our observations. This process re¬ 
vealed evidence that the northern side of the disk exhibits 
a deficit of polarized flux near the inner working angle of 
our data as compared to the southern side of the disk (Fig¬ 
ure ® . We discuss the potential origin of this asymmetry 
in the discussion section. 

4. DISCUSSION 
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Our HiCIAO H-band multi-epoch observations of the 
DoAr 28 transitional disk clearly reveals evidence of a scat¬ 
tered light disk observed in polarized intensity, extending 
from our effective inner working angle of 13 AU (Cf/lO; 2014 
epoch) to 65 AU (0"50). We observe no gap in the disk in 
our PI imagery at the location su ggested by previou s SED 
modeling of the system (15 AU; Kim et al. 2013). This 
suggests that either the small grain population is decou¬ 
pled from the large grain population, and/or that the disk 
gap resides inside of our effective inner working angle. Our 
MCRT modeling of the disk, using constraints both from 
DoAr 28’s SED and our H-band imagery, suggests the disk 
has a gap extending from the dust sublimation radius out 
to ~8 AU, that is only partially cleared of material. This 
gap size is smaller than that esti mated from previo us SED- 
only modeling efforts (15 AU; Kim et al. j2013 ). Other 
transitional disks observed in the H-band with HiCIAO as 
part of the SEEDS project, such as SAP 206462, MWC 


758, and SR 21 (Muto et al.||2012 Grady et al. 

2013 

Fob 

lette et al.|2013 

), also exhibit evidence of small c 

ust grains 


in their gaps, similar to that inferred for DoAr 28. Future 
multi-wavelength observations of the system that achieve 
a factor of ~2 improvement in the effective inner working 
angle are needed to both test if the disk gap is truly as 
small as ~8 AU, and to determine whether the radial dis¬ 
tribution of small and large dust grains are decoupled , as 
is the case for systems like SR 21 (Follette et al.||2013). 

The seemingly small gap size of DoAr 28 contrasts the 
disk gaps at larger orbital separations observed in other 


transitional disks, such as PDS 70 (70 AU; Hashimoto et 
al.|2011l), SR 21 (36 AU;lFollette et al.|2013| ) , SAP 206462 


(46 AU|Muto et al.||2012 ), LkCa 15 (56 AU; Thalmann et 
al.|2014[ ), and Oph IRS 48 (60 AU [Follette et al.||2015|) . R 
this gap is caused by dynamical interactions between plan- 
etary bodies in the system and the disk, this could suggest 
that DoAr 28 has few such fully formed companions, per¬ 
haps owing to the system being in a more youthful state 
than other imaged systems. 

DoAr 28 exhibits several indications of morphological 
features that deviate from simple axisymmetry at modest 
signal-to-noise levels (Figure [^. For example, our 2014 
epoch imagery exhibits tentative evidence of asymmetry 
in the northern-side of the disk (Figure!^. This type of 
asymmetry is commonly observed, and could be caused by 
a range of phenomenon ranging from companion interac¬ 
tions, a clumpy inner disk, and magneto-rotational insta¬ 
bility. A companion such as a planet could induce warping 
of the disk which would shadow the outer portions of the 
disk. If the observed asymmetric feature is caused by a 
companion, the feature should move in the disk at the same 
rotational speed of the companion , and not at the Keplar- 
ian speed ( Hashimoto et al!||2011 ). The asymmetry could 
also indicate the presence of a small azimuthal asymme¬ 
try in the inner disk. Numerical simulations of magneto¬ 
rotations instabilities (MRI) als o suggest that MRI -driven 
disk winds can perturb the disk (Suzuki et al. 2010), which 
could also produce departures from axisymmetry. Addi¬ 
tional epochs of spatially resolved imagery could help con¬ 
strain whether the observed non-axisymmetric structure 
evolves with time, hence constrain the potential origin of 
this feature. 


An additional, albeit more speculative indication of a 
non-axisymmetric feature is present along the southwest 
edge of the disk in our 2014 epoch data, namely the sm- 
gestive presence of a small hook-like feature (Figure [R. 
While the feature looks similar to spiral arms seen on in 


disks such as SAP 206462 and MWC 758 (Muto et al. 


2012 Grady et al. 2013), there are no known T-Tauri tran- 
sitional disks that exhibit spiral structures. It has been 
suggested that the lower temperatures of T-Tauri disks 
(compared to their Herbig analogs) results in lower s ound 
speeds, creating tighter spiral arms ( Muto et al.|20I2 ) that 
would not be resolved in our dataT Some studies have 
found that the dynamical temperature of C O gas is greater 
in T Tauri stars than in Herbig AeBe stars (2013), but this 
CO gas may not be well coupled with the small dust grain 
population. We conclude that this hook, if it is real, is 
most likely to be a perturbation of the disk much like the 
asymmetric feature described above and not a spiral arm. 


5. CONCLUSION 

We have reported the first spatially resolved scattered 
light image of the DoAr 28 transitional disk in H-band. 
We detect the scattered light disk from O'.'IO (13 AU) out 
to 0'.'50 (65 AU), which is slightly interior to the loca¬ 
tion of the system’s gap inferred by another group’s SED 
modeling (15 AU). Although we detected a point source 
companion 171 northwest of the system, our second epoch 
of imagery of the system indicates this object is most likely 
a background star. Using the HOCHUNK3D Monte Carlo 
Radiative Transfer code, we have modeled both the ob¬ 
served SED and H-band PI imagery of the system. Our 
best fit models utilize a modestly inclined (50°), 0.01 M© 
disk that has a partially depleted inner gap from the dust 
sublimation radius out to ~8 AU. Subtracting this best 
fit, axi-symmetric model from our polarized intensity data 
reveals evidence for a small asymmetry in the northern- 
side of the disk, which could be attributable to variety of 
mechanisms. We encourage future high spatial resolution 
sub-millimeter imagery of the system to better ascertain 
the location of the disk gap in the system, and to search 
for azimuthal and radial differences in distribution of small 
versus large dust grains that could be caused by recent 
planet formation in the system. 
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Fig. 1.— The H-band polarized intensity (PI) image of DoAr 28 i n 20 14 June (left) and 2012 July (right) is shown over a field of view of 
F.'95 X l('95 with the excess halo subtracted as described in section [2.1| We have applied a software mask over the location of the central 
star, with a size of CK'17 (24 AU) for the 2012 epoch data and 0^/10 (13 AU) for the 2014 epoch data, for asthetic purposes. The imagery is 
plotted linearly in units of mjy, and was further smoothened with a 3-pixel Gaussian kernel. 


RA Offset (■■) 

0.4 0.2 0 -0.2 -0.4 



Fig. 2.— A polarization vector map is overlaid on the PI image of DoAr 28 of epoch 2014 including halo subtraction. While the direction 
of the vectors represent the orientation of the observed polarization, the length of the vectors is a relative quantity and does not indicate the 
percent polarization present. The clear centrosymmetric behavior of the vectors confirms that we are detecting scattered light originating 
from the system’s disk. 



log(radius) (arcsec) 


Fig. 3.— The radial profile of 2014 epoch data along the north-side of the major axis illustrates clear evidence of two departures from a 
singular power law behavior. One clear break in the power law happens at ~0''l (the FHWM of the 2014 data), where the system transitions 
from being dominated by unsubtracted PSF residuals (blue) to disk dominated scattered light (green). The second break occurs when the disk 
dominated scattered light (green) is overwhelmed by background noise (red). The solid colored lines are best fit lines from linear regression 
fits to the disk power law shown on the figure. The black vertical line represents where the inner working angle is defined. 
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Fig. 4.— The observed surface brightness along the major axis of the disk in our PI imagery at our 2012 (green) and 2014 epochs (blue), 
computed with a 4 pixel-wide average, is shown along with analogous surface brightness measurements from our MCRT model (red). The 
shaded colors represent a 1-sigma error bar with light green for 2012 epoch and light blue for 2014 epoch. The left side corresponds to the 
South end of the disk and the right side corresponds to the North side of the disk shown in Figure^ The central portions of these profiles 
have been removed, signifying the effective inner working angle of our data. 
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Fig. 5.— We subtracted our best fit model from our 2014 epoch PI imagery, yielding the residual imagery shown on the left. The right panel 
represents the signal to noise of the residual imagery. The image has been rotated such that the major axis is horizontal in the figure, with the 
northeast region of the disk on the right-hand side of the image. Two potential asymmetries are observed, including a “hook” feature noted 
in Figure^ and a further asymmetry discussed in SectionThe imagery is plotted linearly in units of ADU, and was further smoothened 
with a 3-pixel Gaussian kernel. 
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Fig. 6. — The proper motion of DoAr 28 (blue star) and its candidate companion (red circle) in our two epochs of observations. DoAr 
28 shares proper motion with p Oph shown as the b lack sinusoidal line. The observed proper motion suggests that the companion is not 
CO- located with DoAr 28, as discussed in Section[3.2| 


60 AU 100 AU 200 AU 300 AU 



Projected Separation 


Fig. 7. — The computed sensitivity limits of point source detections in our DoAr 28 imagery with a signal to noise of 5 sigma are given, 
along with associated planetary mass limits assuming both BT-Settl (burgundy) and SB12 (green) planet formation models. We assume an 
age of 5 Myr and SB12 “warm start” models midway between the minimum and maximum available initial entropies. We do not detect 
evidence of any co-moving companions above these limits. 
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Fig. 8. — The observed SED of DoAr 28, as compiled in Tablef^ is shown along with our best fit model, using the parameters compiled in 
Table The solid line represents the total SED, whereas the dashed lines represent the different contributions to the SED, including thermal 
emission from the dust disk (thermal), the stellar photosphere (stellar), and scattered light from the disk (scattered). 


Table 1 

Disk Model Parameters 


Variable 

Model Low Bound Upper Bound 

A// (K) 

4375 



Rq 

1.6 

1.5 

1.6 

Disk Mass (Mq) 

0.01 

0.005 

0.05 

Max Disk Radius (AU) 

250 

200 

500 

Fraction of Mass 

0.8 

0.7 

0.99 

Inner Gap Radius (AU) 

0 



Outer Gap Radius (AU) 

8 

7 

9 

Wall height (AU) 

0.5 

0.5 

0.5 

Wall length (AU) 

0.5 

0.5 

1.1 

Zscale Disk 1 

1.7 



Zscale Disk 2 

1.7 

1.7 

1.9 

a 

2.0 

1.9 

2.1 

p 

1.0 

0.99 

1.01 

Accretion {Mq) 

4.0E-9 



Gap Density 

5E-6 

5E-5 

5E-7 

Inclination (i) 

50° 

40° 

65° 


Note. — A summary of the lower and upper acceptable bounds of key HOCHUNK3D parameters that yield a fit consistent with our 
observations, along with the “best” values that were adopted for our final model. The “fraction of mass” parameter represents the fraction 
of the mass located in the large dust grain disk; the rest of the disk mass is located in the small dust grain disk. The parameter “zscale” 
represents the scale height parameter. The “gap density” is the ratio of the dust density inside the gap relative to the density of dust at the 
inner edge of the disk. The ” Zscale Disk” is the scale hight of disks 1 and 2 where disk 1 is large grain dust and disk 2 is small grain dust. 
Parameters ol and describe the density profile of the disk defined in equation^ 












Table 2 

Tabulated SED 


Wavelength (/tm) 

Flux (Jy) 

Error (Jy) 

Facility 


Reference 

0.44 

0.0098 


Mt. Maidanak Obs. ( 

Grankin et al. 

2007 

I 

0.55 

0.0112 


Mt. Maidanak Obs. ( 

Grankin et al. 

20UV 

I 

0.64 

0.0411 


USNO 


Vlonet et al. 

1998 

I 

1.26 

0.1759 

0.0039 

2MASS 


i Gutri et al. 

2003 

I 

1.6 

0.2606 

0.0055 

2MASS 


(Gutri et al. 

2003 

I 

2.22 

0.2379 

0.0042 

2MASS 


(Gutri et al. 

2003 

I 

3.4 

0.1670 

0.0015 

WISE 

Gutri & et al. 

2012 

I 

4.6 

0.1092 

0.00087 

WISE 

Gutri & et al. 

2012 

I 

12.0 

0.0471 

0.00045 

WISE 

Gutri & et al. 

2012 

I 

22.0 

0.4450 

0.0037 

WISE 

Gutri & et al. 

2012 

I 

24.0 

0.348 

0.0322 

Spitzer/MIPS 

(lEvans et al. 

2003 

I 

62 

0.829 

0.013 

Herschel 

( 

Eeane et al. 

2014 

I 

70.0 

0.579 

0.0891 

Spitzer/MIPS 

( 

Evans et al. 

2003 

I 

1300 

0.072 

0.0019 

SMA 


SMA Observation 


The tabulated version of DoAr 28’s photometry that was used to construct the SED shown in Figure]^ 






































